Stability of xenon oxides at high pressures 
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Xenon, which is quite inert under ambient conditions, may become reactive under pressure. The 
possibility of formation of stable xenon oxides and silicates in the interior of the Earth could explain 
the atomspheric missing xenon paradox. Using the ab initio evolutionary algorithm, we predict 
the thermodynamical stabilization of Xe-0 compounds at high pressures (XeO, Xe02 and XeOs at 
pressures above 83, 102 and 114 GPa, respectively). Our calculations indicate large charge transfer 
in these oxides, suggesting that large electronegativity difference and pressure are the key factors 
favoring the formation of xenon compounds. Xenon compounds in the Earth's mantle, however, 
cannot directly explain the missing xenon paradox: xenon oxides are unstable in equilibrium with 
metallic iron in the Earth's lower mantle, while xenon silicates are predicted to spontaneously 
decompose at all mantle pressures (<136 GPa). This does not preclude Xe atoms from being 
retained in defects of mantle silicates and oxides. 



Xenon is a noble gas, chemically inert at ambient con- 
ditions. A few xenon fluorides have been found [IH3], 
with Xe atoms in the oxidation states +2, +4, +6, or 
+8. Upon application of high pressure, the molecular 
phase of insulating XeF2 has been reported to transform 
into two- and three-dimensional extended solids and be- 
come metallic [1]. Clathrate Xe-H solids were also ob- 
served at pressures above 4.8 GPa [5 . Two xenon oxides 
(Xe03, Xe04) [6] are known at atmospheric pressure, 
but are unstable and decompose explosively above 25° 
C (Xe0 3 ) and -40° C (Xe0 4 ) [7 . A crystalline Xe0 2 
phase with local square-planar Xe04 geometry has re- 
cently been synthsized at ambient conditions [5]. 

There is growing evidence to suggest that noble gases, 
especially Xe, may become reactive under pressure [?J]. 
The formation of stable xenon oxides and silicates could 
explain the missing xenon paradox - the observation that 
the amount of Xe in the Earth atmosphere is an order 
of magnitude less than what it would be if all Xe were 
degassed from the mantle into the atomosphere |10j.One 
explaination for this deficit is that Xe is largely retained 
in the mantle. In fact, it has been reported that a few 
weight percent of xenon can be incorporated into Si02 
at elevated pressures and high temperatures pTj H2] . A 
recent theoretical investigation has shown that no xenon 
carbides are stable up to at least 200 GPa [T3] , and an ex- 
perimental and theoretical high pressure study [14] found 
no tendency for xenon to form a metal alloy with iron or 
platinum. 

Here we investigate the possible stability of xenon ox- 
ides using quantum-mechanical calculations of their en- 
ergetics. As the structures of stable xenon oxides are not 
experimentally known, we calculate them using the re- 



cently developed evolutionary algorithm for crystal struc- 
ture prediction [131 HI] • We also analyse chemical bond- 
ing in these exotic compounds. 

We have performed structure prediction simulations 
for the Xe-0 system for the compositions of XeO, Xe02, 
Xe0 3 , Xe0 4 at 5, 50, 100, 120, 150, 180, 200 and 220 
GPa. Our calculation at 5 GPa yielded lowest-enthalpy 
structures that always contained the O2 molecules, indi- 
cating the tendency for segregation of the elements, and 
indeed at 5 GPa decomposition was found to be ener- 
getically favourable. This suggests that the reaction ob- 
served by Sanloup et al. (TTJH2] at 0.7 to 10 GPa was an 
entropically driven incorporation of Xe impurities into 
the structure of Si02, rather than enthalpically-driven 
formation of a stoichiometric xenon silicate or oxide. In- 
deed, solid solutions and point defects are stabilized by 
entropy [T7] . 

However, detailed calculations at and above 100 GPa 
indeed found a number of oxides stable against decom- 
position into the elements. The most promising struc- 
tures and compositions were investigated further. The 
Gibbs free energy of formation of these oxides at high P, 
T conditions is shown in Fig. [T] At K, XeO, Xe0 2 
and Xe03 at high P become thermodynamically stable. 
We also performed quasi-harmonic free energy calcula- 
tions with phonon spectra using the finite-displacement 
method [18] to assess the effect of tcmperture on their 
stability (Fig. ^p) . The inclusion of thermal effects does 
not bring any quantitative changes into the picture, and 
XeO, Xe02 and Xe03 remain thermodynamically stable 
also at high T. 

Fig. [2]i shows the enthalpy of formation of all the 
Xe oxides as a function of pressure. Below 83 GPa all 
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FIG. 1. a) Predicted enthalpy of formation of Xe-O compounds at high P and T=Q K; b) Predicted Gibbs free energy of 
formation of Xe-O compounds at different temperatures and P=120 GPa. The compounds shown are, from left to right, Xe, 
XeO, Xe02, Xe03, O. For oxygen, the structures of the ("-phase |19] and e-phase [20] were used. For Xe, the fee |21j and hep 
|22| structures were considered and hep was found energetically more favourable above 100 GPa, in agreement with experiments 



xenon oxides are unstable. At 83 GPa, XeO-Pbcm be- 
comes stable, followed by Xe02-P2i/c above 102 GPa 
and Xe03-P42/mnm above 114 GPa. There is an inter- 
esting trend for the oxidation number of Xe to increase 
with increasing pressure. 

A simple and clear analysis of chemical bonding can 
be done using the electron localization function (ELF) 
[21]. The ELF gives information about the valence elec- 
tron configuration of an atom in a compound. States 
with closed shell electron configurations(Xe , 5s 2 5p 6 , and 
Xe 6+ , 5s 2 ) will exhibit a spherical ELF distribution, 
whereas open shell states (Xe 2+ , Xe 4+ ) will not. For 
Xe 2+ one p-orbital is empty and the ELF will have a 
toroidal shape; likewise, Xe 4+ can be formed by the re- 
moval of two p-orbitals and the ELF will show a two- 
lobe maximum corresponding to the shape of the lone 
p-electron pair. 

The most stable structure of XeO at 100 GPa has space 
group Pbcm and 8 atoms in the unit cell. As shown in 
Fig. [2^,, Xe atoms are in a twofold (linear) coordina- 
tion and Xe-O bonds form chains, with O-Xe-0 angles 
of 175.6° and Xe-O-Xe angles of 112.6°. The alternating 
Xe-O bond lengths are 2.0 and 2.1 A. The ELF picture 
(Fig. [2^,) shows a toroidal maximum of ELF around each 
Xe atom, exactly what one should expect for Xe 2+ state. 
Above 145 GPa XeO undergoes a phase transformation 
and forms a structure with the space group P2\/ m. Xel 
connects to four oxygens and has square coordination, 
and forms thus the same chains as in Xe02 (suggest- 
ing that Xel atoms are in the tetravalent Xe 4+ state), 
whereas Xe2 can be described as neutral and not bonded 
to other atoms by any significant bonds. The presence 
of neutral non-bonded atoms in this structure is energet- 
ically favourable as it increases its packing density. 

For Xe02, the stable structure above 102 GPa has 
space group P2\j c and 24 atoms in the unit cell. Xenon 
atoms have a slightly non-planar square coordination 



and the structure consists of ID-ribbons of edge-sharing 
Xe04-squares (Xe-O distances are 2.0 and 2.1 A), with 
four Xe-O bonds and two lone-pair maxima forming an 
octahedron, consistent with the geometry proposed by re- 
cent experiment [§]. Just as in XeO, there are no peaks 
visible in the ELF isosurface along the Xe-O bonds (Fig. 
[2jj). Above 198 GPa Xe0 2 transforms into the Xe0 2 - 
Cmcm structure. This non-trivial structure can be rep- 
resented as having parallel zigzag Xe chains (Xe-Xe dis- 
tances are 2.62 A), with each Xe atom having two neigh- 
bouring oxygen atoms in the form reminiscent of bent 
(the Xe-O distance is 1.95 A, O-Xe-0 angle is 160.7°) 
Xe02 molecules. 

Xe03 becomes stable at 114 GPa. Its structure has 
space group is P^/mnm and 16 atoms in the unit cell. 
It is stable towards decomposition into Xe and O2 as well 
as into XeO or Xe02 and O2. As shown in Supplemen- 
tary Materials, P^/mnm phase is composed of two sub- 
lattices: square Xe02 chains again, suggesting the Xe 4+ 
state and chains made of linear O2 dumbbells. Above 
145 GPa, the molecules in the linear -O2-O2- chains are 
partly dissociated and we observe the -O2-O- chain in the 
C2/ c phase with 48 atoms per unit cell. Above 198 GPa, 
the structure transforms to a Pmmn phase with 8 atoms 
per unit cell (Supplementary Materials). In this remark- 
able structure, the oxygen atoms form anticuboctahedrea 
in which the Xe atom sits in the center (Fig. [2J2). The 
ELF distribution around Xe atoms in the Pmmn phase is 
spherical around the xenon, which points at the Xe 6+ va- 
lence state with a spherically symmetric 5s 2 valence shell. 
Again, we observe the tendency of increasing oxidation 
states under pressure. 

Xenon fluorides are stable at ambient conditions, 
xenon oxides become stable above 83 GPa, and xenon 
carbides are unstable up to 200 GPa at least [13] ■ It thus 
appears that xenon forms compounds most readily with 
the most electronegative atoms, and that in turn suggests 



(a) 



(c) 



» - • • • 






XeO-Pbcm 

- - XeO-Pl^m 

— Xe0 2 -P^-[lc 
XeOyPtymnm 
Xe0 3 -C2/c 



FIG. 2. a) Crystal structure of XeO (Pbcm) Structure at 100 GPa, and its ELF picture (p = 0.85) on the Xe-O chain; b) 
Crystal structure of Xe02 (P2\/c) structure at 120 GPa, and its ELF picture (p — 0.85) on the Xe04 square, c) crystal 
structure of Xe03 (Pmmn) Structure at 200 GPa, and its ELF picture (p = 0.82) on XeOi2 anticuboctahedra; d) Enthalpies 
of formation (solid lines) and enthalpies of decomposition into other xenon oxides (dashed lines) of all stable xenon oxides. 
Stability fields are delineated by colors (At 198 GPa, P2\j c Xe02 transforms to a Cmcm structure, and C2/c-XeOa transforms 
to Pmmn phase, i.e., just at the edge of the graph. 



that ionicity is essential. This is somewhat counterintu- 
itive, given that xenon atom has a very stable closed va- 
lence shell and its electronegativity is rather high. The 
electronegativity difference (1.4 for Xe-F, 0.8 for Xe-0 
and 0.56 for Xe-C) determines the degree of ionicity at 
ambient conditions. However, ionicity could be enhanced 
under pressure. Spontaneous ionization under pressure 
was recently found even in elemental boron |25] . 

To obtain further chemical insight into these exotic 
xenon oxides, we selected some of the stable structures 
containing Xe 2+ , Xe 4+ , and Xe 6+ , namely, XeO- Pbcm 
at 100 GPa, Xe0 2 -P2i/c at 150 GPa, and Xe0 3 Pmmn 
at 200 GPa. Fig. [3] shows the density of states and 
its projection onto atomic orbitals. All of these xenon 
oxides are narrow-gap semiconductors. Using state-of- 
the-art GW calculations, we got the band gaps: 1.52 
eV for XeO-Pbcm, 0.52 eV for Xe0 2 -P2i/c, 0.15 eV for 
XeO^,-Pmmn; these values should be accurate to within 
5-10%. 

The highest valence band levels are dominated by p 
orbitals of O and Xe. Both the atom-projected densities 
of states and energy-decomposed electron densities sug- 
gest charge transfer from Xe to O atoms. Note that the 
contribution of Xe p orbitals strongly decreases with in- 
creasing Xe oxidation because of two concurring effects, 



the change in stoichiometry and the parallel enhanced 
electron transfer from Xe to O. These two effects are 
disentangled in the DOS reported in Supplementary Ma- 
terial. The valence states dominated by p contributions 
contain in XeO-Pbcm about 40 electrons, with almost 
equal contributions from O and Xe. The contribution is 
almost all p, namely 4.51 p electron per Xe and 4.93 p 
electrons per O atom. In the case of Xe02--P2i/c there 
are about 112 electrons, 31.94 from Xe and 80.16 from 
O, which implies about 4 e per Xe and 5 e per O, which 
suggests that about one electron is transferred from a 
p-orbital of Xe into a p-orbital of O. Finally, in Xe03 
Pmmn, there are about 45 electrons in the valence states 
dominated by p orbital contributions, 5.48 e from Xe and 
30.10 e from O, leading to a further lowering to 2.8 elec- 
trons on Xe p orbitals and to again about 5 electrons on 
O p orbitals. 

The charge transfer issue was also investigated on the 
basis of the electron density using Bader's analysis |26j . 
For XeO-Pbcm the net charge on the Xe atoms is 1.01, 
while it increases to 1.995 on average in Xe02-P2i/c and 
to 2.754 in XeO^-Pmrnn. Note that the net charge on 
oxygen atoms stay almost the same in all the three com- 
pounds and very close to -1. This corresponds to ionicity 
of about 50% in all oxides. Note that the charge transfer 
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FIG. 3. Total and projected densities of states for a) XeO-P&cm at 100 GPa, b) Xe02--P2i/c at 150 GPa, c) 'KeOs-Pmmn 
at 200 GPa. The sum over the plotted projections gives the total DOS for each system. The DOS plotted have different units 
on the y-axis so as to give the total number of valence electrons for each system when integrated up to the Fermi level. This 
corresponds to 56 e" (XeO, Z=4), 160 e" (Xe0 2 , Z=8) and 52 e" (Xe0 3 , Z=2). 



described through the observable function, the electron 
density, and Bader's analysis confirms what was found 
with using orbital-projected DOS, a totally independent 
approach. 

Bader analysis yields not only the net atomic charges, 
but also more subtle characteristics such as the three 
eigenvalues /ij (i=l,3) of the traceless quadrupole mo- 
ment tensor Q(fi!), which is a good quantitative indica- 
tor of the departure of an atomic basin from sphericity. 
The eigenvectors associated with the eigenvalues /ij give 
the principal directions for relative charge accumulation 
and depletion. For a spherical distribution the fj,^ are all 
equal to zero. Deviation from zero indicates asphericity: 
negative eigenvalues /ij arise from accumulation of charge 
in the direction associated with the corresponding eigen- 
vector and at the expense of the directions associated 
with positive fij. Our calculation (Table]!]) shows that Xe 
atoms in HsO^-Pmmn are definitely more spherical, their 
/ij being all close to zero and about one order of magni- 
tude lower than for the Xe atoms in Xe02 and XeO. This 
agrees with the ELF picture while providing also a quan- 
titative measure of sphericity. XeO-P&cm (with one p 
orbital partially empty) has, in agreement with the or- 
bital picture, one direction of relative charge depletion, 
associated to /13, and two directions of unequal relative 
charge accumulation, associated to fi\ and \xi- Finally in 
Xe02--P2i/c, there is one direction of relative charge ac- 
cumulation associated to fj,i and two of unequal relative 
charge depletion, associated to 112 and /Z3, in agreement 
with the picture of two p orbitals (partially) empty. 

How likely is it that xenon oxides (or silicates) could 
exist in the Earth's lower mantle, thus explaining the 
missing Xe paradox? As discussed above, Xe oxides are 
only stable above 83 GPa, i.e. at pressures corresponding 



to the lower mantle, where metallic Fe should be present 
[27l [28] . All of the predicted xenon oxides are very strong 
oxidants and our calculations show that they will be re- 
duced by Fe, producing iron oxide and free xenon. Then 
we investigated the formation of xenon silicates under 
pressure, focusing on XeSiOa and Xe2Si04, which con- 
tain the least oxidized xenon. All of the investigated 
compositions were unstable towards decomposition into 
XeO, Xe02, Si02 and the elemental Xe; Xe2Si04 proved 
to be one of the least unstable silicates, but still is un- 
stable. Our conclusion is: while Xe oxides are stable in 
the Xe-0 system, at strongly reducing conditions neither 
silicates nor oxides of Xe can exist: Xe silicates are un- 
stable to decomposition, while Xe oxides will be reduced 
by metallic Fe present in the lower mantle. 

To summarize, we have predicted stability of xenon 
oxides at high pressure, which can be readily tested ex- 
perimentally. On increasing pressure, increasingly high 
oxidation states of Xe will appear - first, XeO (above 83 
GPa), then Xe0 2 (above 102 GPa), then Xe0 3 (above 
114 GPa). Present results clearly show that Xe loses 
its chemical inertness under pressure and that charge 
transfer plays an essential role in chemical bonding of 
Xe compounds and their stability is largely determined 
by electronegativity differences. Furthermore, pressure 
stabilizes increasing oxidation states of Xe atoms (Xe° 
— ► Xe 2+ — ► Xe 4+ — > Xc 6+ ) and enhances charge trans- 
fer from Xe to O atoms. We find that xenon silicates are 
not stable at pressures of the Earth's mantle (<136 GPa) 
and that xenon oxides, although stable against decompo- 
sition into the elements, will be reduced to free Xe in the 
strongly reducing conditions of the lower mantle. 

Although the formation of stable Xe oxides or sili- 
cates is not possible at conditions found in the Earth's 
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System 


Mi 


M2 ^3 


XeO 


-3.49 


-2.04 5.53 


Xe0 2 -Xel 


-5.50 


1.83 3.67 


Xe0 2 -Xe2 


-5.90 


2.33 3.57 


XeQ 3 


-0.38 


0.21 0.36 



TABLE I. Eigenvalues of the traceless quadrupole moment tensor from the atomic Bader intergration. The selected xenon 
oxide structures are: XeO-Pbcm at 100 GPa, Xe0 2 -P2i/c at 150 GPa, and Xe0 3 -Pmmn at 200 GPa 



mantle, the formation of strong Xe-0 bonds under pres- 
sure, clearly seen in our results, implies that Xe may 
still be retained at point or line defects, or grain bound- 
aries of mantle minerals. Xenon could also be stored in 
perovskite/post-perovskite stacking faults [33]. The facile 
chemical bonding between Xe and O atoms demonstrated 
here and the preference of Xe atoms to terminate the sili- 
cate perovskite layers, observed in our simulations (xenon 
silicate in Supplementary Materials), both suggest this 
possibility. Indeed, the effect of trapping of trace ele- 
ments by lattice defects is well known [3U] 

METHODS 

Crystal Structure Prediction: The evolutionary 
algorithm USPEX [HI [16] , used here for predicting new 
stable structures, searches for the structure with the low- 
est free energy at the given pressure-temperature condi- 
tions and is capable of predicting the stable structure 
of a compound knowing just the chemical composition. 
Details of the method are described elsewhere [El US] 
and a number of applications [HI [25] illustrate its power. 
All structure prediction runs discussed here were done 
with the USPEX code, with structure relaxations done 
using density functional theory (DFT) within the gener- 
alized gradient approximation (GGA) [31] in the frame- 
work of all-electron projector augmented wave (PAW) 
[32"] method as implemented in the VASP [Sj^ code. We 
used the plane wave kinetic energy cutoff of 520 eV and 
the Brillouin zone was sampled with the resolution of 27r 
x 0.08 A -1 , which showed excellent convergences of the 
energy differences, stress tensors and structural param- 
eters. We studied systems containing up to 36 atoms 
per unit cell. The first generation of structures was cre- 
ated randomly. All structures were relaxed at constant 
pressure and the enthalpy was used as fitness. The en- 
ergetically worst structures (40%) were discarded and a 
new generation was created from the remaining struc- 
tures through heredity, lattice mutation and permutation 
of atoms. Additionally, the best structure of a genera- 
tion was carried over into the next generation. We ter- 
minated the runs generally after 50 generations, and all 
runs had found the minimum enthalpy structures much 
earlier. The population size was set to at least twice 
the number of atoms in the cell. Results obtained with 
and without the van der Waals functional 34 (as imple- 



mented in the VASP code) are very similar, and here we 
show the results that include this functional. 

Chemical bonding analysis: The density of states 
and its projection onto the atomic orbitals have been cal- 
culated through the periodic LCAO approach, using the 
CRYSTAL-06 code [35] and the same DFT functional [3T] 
as used in all calculations described here. Oxygen atoms 
were described by an 8-411+d all electron basis set [3J5] 
and Xe atoms by a cc-pVTZ basis set [37], specifically 
devised for the fully-relativistic ECP28MDF pseudopo- 
tential, with 26 electrons per Xe kept active. In order to 
avoid problems of numerical catastrophes |35j , the origi- 
nal cc-pVTZ(12sllp9dl/)/[5s4p3(il/] basis set [37] was 
slightly modified by removing the outermost s, p and d 
Gaussians and the / polarization function to yield a final 
(Ilsl0p8d)/[5s4p3d] contracted basis set. The Kohn- 
Sham matrix was diagonalized on an isotropic 8x8x8 
k-mesh and the same mesh was used in the Fermi energy 
calculation and density matrix reconstruction. Bader's 
charges and atomic quadrupole moment tensors were 
evaluated from the periodic LCAO electron density and 
using the TOPOND package [35] interfaced to the CRYS- 
TAL code. TOPOND determines the boundaries of the 
atomic basins and the integrated properties within these 
basins using fully analytical and on-the-fiy evaluations 
of the electron densities and its derivatives (up to the 
4th order), i.e. no use is made of electron densities on a 
grid or of numerical approximations of electron density 
gradients |39j . Basin boundaries are determined through 
the PROMEGA algorithm @D], while the basin integra- 
tion is performed in spherical coordinates, using Gaus- 
sian quadrature formulas. The resulting atomic charges 
were very close to those obtained from the PAW wave- 
function and the numerical algorithm described in [41) . 
Band gaps were recalculated within GW approximation 
as implemented in the VASP code [32"] . 

Calculations were performed on the CFN cluster and 
Blue Gene supercomputer (Brookhavcn National Labo- 
ratory), Swiss Supercomputer Centre, Skif MSU super- 
computer (Moscow State University) and at the Joint 
Supercomputer Center of the Russian Academy of Sci- 
ences (Moscow). A.R.O. thanks DARPA (grant no. 
W31P4Q1210008) and the National Science Foundation 
(grant no. EAR-1114313) for financial support. 
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